Normal human mammary epithelial cells spontaneously escape senescence and acquire genomic changes

Senescence and genomic integrity are thought to be important barriers in the development of malignant lesions. Human fibroblasts undergo a limited number of cell divisions before entering an irreversible arrest, called senescence. Here we show that human mammary epithelial cells (HMECs) do not conform to this paradigm of senescence. In contrast to fibroblasts, HMECs exhibit an initial growth phase that is followed by a transient growth plateau (termed selection or M0; refs 3-5), from which proliferative cells emerge to undergo further population doublings (approximately 20-70), before entering a second growth plateau (previously termed senescence or M1; refs 4-6). We find that the first growth plateau exhibits characteristics of senescence but is not an insurmountable barrier to further growth. HMECs emerge from senescence, exhibit eroding telomeric sequences and ultimately enter telomere-based crisis to generate the types of chromosomal abnormalities seen in the earliest lesions of breast cancer. Growth past senescent barriers may be a pivotal event in the earliest steps of carcinogenesis, providing many genetic changes that predicate oncogenic evolution. The differences between epithelial cells and fibroblasts provide new insights into the mechanistic basis of neoplastic transformation.

Integrated analysis of epigenomic and genomic changes by DNA methylation dependent mechanisms provides potential novel biomarkers for prostate cancer

Epigenetic silencing mediated by CpG methylation is a common feature of many cancers. Characterizing aberrant DNA methylation changes associated with tumor progression may identify potential prognostic markers for prostate cancer (PCa). We treated two PCa cell lines, 22Rv1 and DU-145 with the demethylating agent 5-Aza 2’–deoxycitidine (DAC) and global methylation status was analyzed by performing methylation-sensitive restriction enzyme based differential methylation hybridization strategy followed by genome-wide CpG methylation array profiling. In addition, we examined gene expression changes using a custom microarray. Gene Set Enrichment Analysis (GSEA) identified the most significantly dysregulated pathways. In addition, we assessed methylation status of candidate genes that showed reduced CpG methylation and increased gene expression after DAC treatment, in Gleason score (GS) 8 vs. GS6 patients using three independent cohorts of patients; the publically available The Cancer Genome Atlas (TCGA) dataset, and two separate patient cohorts. Our analysis, by integrating methylation and gene expression in PCa cell lines, combined with patient tumor data, identified novel potential biomarkers for PCa patients. These markers may help elucidate the pathogenesis of PCa and represent potential prognostic markers for PCa patients.

Rare genomic structural variants in complex disease: lessons from the replication of associations with obesity.

The limited ability of common variants to account for the genetic contribution to complex disease has prompted searches for rare variants of large effect, to partly explain the 'missing heritability'. Analyses of genome-wide genotyping data have identified genomic structural variants (GSVs) as a source of such rare causal variants. Recent studies have reported multiple GSV loci associated with risk of obesity. We attempted to replicate these associations by similar analysis of two familial-obesity case-control cohorts and a population cohort, and detected GSVs at 11 out of 18 loci, at frequencies similar to those previously reported. Based on their reported frequencies and effect sizes (OR≥25), we had sufficient statistical power to detect the large majority (80%) of genuine associations at these loci. However, only one obesity association was replicated. Deletion of a 220 kb region on chromosome 16p11.2 has a carrier population frequency of 2×10(-4) (95% confidence interval [9.6×10(-5)-3.1×10(-4)]); accounts overall for 0.5% [0.19%-0.82%] of severe childhood obesity cases (P = 3.8×10(-10); odds ratio = 25.0 [9.9-60.6]); and results in a mean body mass index (BMI) increase of 5.8 kg.m(-2) [1.8-10.3] in adults from the general population. We also attempted replication using BMI as a quantitative trait in our population cohort; associations with BMI at or near nominal significance were detected at two further loci near KIF2B and within FOXP2, but these did not survive correction for multiple testing. These findings emphasise several issues of importance when conducting rare GSV association, including the need for careful cohort selection and replication strategy, accurate GSV identification, and appropriate correction for multiple testing and/or control of false discovery rate. Moreover, they highlight the potential difficulty in replicating rare CNV associations across different populations. Nevertheless, we show that such studies are potentially valuable for the identification of variants making an appreciable contribution to complex disease.

In planta conditions induce genomic changes in Pseudomonas syringae pv. phaseolicola

The co-evolution of bacterial plant pathogens and their hosts is a complex and dynamic process. Plant resistance can impose stress on invading pathogens that can lead to, and select for, beneficial changes in the bacterial genome. The Pseudomonas syringae pv. phaseolicola (Pph) genomic island PPHGI-1 carries an effector gene, avrPphB (hopAR1), which triggers the hypersensitive reaction in bean plants carrying the R3 resistance gene. Interaction between avrPphB and R3 generates an antimicrobial environment within the plant, resulting in the excision of PPHGI-1 and its loss from the genome. The loss of PPHGI-1 leads to the generation of a Pph strain able to cause disease in the plant. In this study, we observed that lower bacterial densities inoculated into resistant bean (Phaseolus vulgaris) plants resulted in quicker PPHGI-1 loss from the population, and that loss of the island was strongly influenced by the type of plant resistance encountered by the bacteria. In addition, we found that a number of changes occurred in the bacterial genome during growth in the plant, whether or not PPHGI-1 was lost. We also present evidence that the circular PPHGI-1 episome is able to replicate autonomously when excised from the genome. These results shed more light onto the plasticity of the bacterial genome as it is influenced by in planta conditions.

Toward resolving deep neoaves phylogeny : data, signal enhancement, and priors

We report three developments toward resolving the challenge of the apparent basal polytomy of neoavian birds. First, we describe improved conditional down-weighting techniques to reduce noise relative to signal for deeper divergences and find increased agreement between data sets. Second, we present formulae for calculating the probabilities of finding predefined groupings in the optimal tree. Finally, we report a significant increase in data: nine new mitochondrial (mt) genomes (the dollarbird, New Zealand kingfisher, great potoo, Australian owlet-nightjar, white-tailed trogon, barn owl, a roadrunner [a ground cuckoo], New Zealand long-tailed cuckoo, and the peach-faced lovebird) and together they provide data for each of the six main groups of Neoaves proposed by Cracraft J (2001). We use his six main groups of modern birds as priors for evaluation of results. These include passerines, cuckoos, parrots, and three other groups termed “WoodKing” (woodpeckers/rollers/kingfishers), “SCA” (owls/potoos/owlet-nightjars/hummingbirds/swifts), and “Conglomerati.” In general, the support is highly significant with just two exceptions, the owls move from the “SCA” group to the raptors, particularly accipitrids (buzzards/eagles) and the osprey, and the shorebirds may be an independent group from the rest of the “Conglomerati”. Molecular dating mt genomes support a major diversification of at least 12 neoavian lineages in the Late Cretaceous. Our results form a basis for further testing with both nuclear-coding sequences and rare genomic changes.

A mitochondrial genome phylogeny of termites (Blattodea: Termitoidae) : robust support for interfamilial relationships and molecular synapomorphies define major clades

Despite their ecological significance as decomposers and their evolutionary significance as the most speciose eusocial insect group outside the Hymenoptera, termite (Blattodea: Termitoidae or Isoptera) evolutionary relationships have yet to be well resolved. Previous morphological and molecular analyses strongly conflict at the family level and are marked by poor support for backbone nodes. A mitochondrial (mt) genome phylogeny of termites was produced to test relationships between the recognised termite families, improve nodal support and test the phylogenetic utility of rare genomic changes found in the termite mt genome. Complete mt genomes were sequenced for 7 of the 9 extant termite families with additional representatives of each of the two most speciose families Rhinotermitidae (3 of 7 subfamilies) and Termitidae (3 of 8 subfamilies). The mt genome of the well supported sister group of termites, the subsocial cockroach Cryptocercus, was also sequenced. A highly supported tree of termite relationships was produced by all analytical methods and data treatment approaches, however the relationship of the termites + Cryptocercus clade to other cockroach lineages was highly affected by the strong nucleotide compositional bias found in termites relative to other dictyopterans. The phylogeny supports previously proposed suprafamilial termite lineages, the Euisoptera and Neoisoptera, a later derived Kalotermitidae as sister group of the Neoisoptera and a monophyletic clade of dampwood (Stolotermitidae, Archotermopsidae) and harvester termites (Hodotermitidae). In contrast to previous termite phylogenetic studies, nodal supports were very high for family-level relationships within termites. Two rare genomic changes in the mt genome control region were found to be molecular synapomorphies for major clades. An elongated stem-loop structure defined the clade Polyphagidae + (Cryptocercus + termites), and a further series of compensatory base changes in this stem loop is synapomorphic for the Neoisoptera. The complicated repeat structures first identified in Reticulitermes, composed of short (A-type) and long (B-type repeats) defines the clade Heterotermitinae + Termitidae, while the secondary loss of A-type repeats is synapomorphic for the non-macrotermitine Termitidae.

Cross-species chromosome painting in the golden mole and elephant-shrew: support for the mammalian clades Afrotheria and Afroinsectiphillia but not Afroinsectivora

Cross-species painting (fluorescence in situ hybridization) with 23 human (Homo sapiens (HSA)) chromosome-specific painting probes (HSA 1-22 and the X) was used to delimit regions of homology on the chromosomes of the golden mole (Ghrysochloris asiaticus) and elephant-shrew (Elephantulus rupestris). A cladistic interpretation of our data provides evidence of two unique associations, HSA 1/19p and 5/21/3, that support Afrotheria. The recognition of HSA 5/3/21 expands on the 3/21 synteny originally designated as an ancestral state for all eutherians. We have identified one adjacent segment combination (HSA2/8p/4) that is supportive of Afroinsectiphillia (aardvark, golden mole, elephant-shrew). Two segmental combinations (HSA 10q/17 and HSA 3/20) unite the aardvark and elephant-shrews as sister taxa. The finding that segmental syntenies in evolutionarily distant taxa can improve phylogenetic resolution suggests that they may be useful for testing sequence-based phylogenies of the early eutherian mammals. They may even suggest clades that sequence trees are not recovering with any consistency and thus encourage the search for additional rare genomic changes among afrotheres.

Comparative cytogenetics of bats (Chiroptera): The prevalence of Robertsonian translocations limits the power of chromosomal characters in resolving interfamily phylogenetic relationships

Although the monophyly of Chiroptera is well supported by many independent studies, higher-level systematics, e.g. the monophyly of microbats, remains disputed by morphological and molecular studies. Chromosomal rearrangements, as one type of rare genomic changes, have become increasingly popular in phylogenetic studies as alternatives to molecular and other morphological characters. Here, the representatives of families Megadermatidae and Emballonuridae are studied by comparative chromosome painting for the first time. The results have been integrated into published comparative maps, providing an opportunity to assess genome-wide chromosomal homologies between the representatives of eight bat families. Our results further substantiate the wide occurrence of Robertsonian translocations in bats, with the possible involvement of whole-arm reciprocal translocations (WARTs). In order to search for valid cytogenetic signature(s) for each family and superfamily, evolutionary chromosomal rearrangements identified by chromosomal painting and/or banding comparison are subjected to two independent analyses: (1) a cladistic analysis using parsimony and (2) the mapping of these chromosomal changes onto the molecularly defined phylogenetic tree available fromthe literature. Both analyses clearly indicate the prevalence of homoplasic events that reduce the reliability of chromosomal characters for resolving interfamily relationships in bats.

Genomic and functional profiling of gastric cancer

Helicobacter pylori infection is a risk factor for gastric cancer, which is a major health issue worldwide. Gastric cancer has a poor prognosis due to the unnoticeable progression of the disease and surgery is the only available treatment in gastric cancer. Therefore, gastric cancer patients would greatly benefit from identifying biomarker genes that would improve diagnostic and prognostic prediction and provide targets for molecular therapies. DNA copy number amplifications are the hallmarks of cancers in various anatomical locations. Mechanisms of amplification predict that DNA double-strand breaks occur at the margins of the amplified region. The first objective of this thesis was to identify the genes that were differentially expressed in H. pylori infection as well as the transcription factors and signal transduction pathways that were associated with the gene expression changes. The second objective was to identify putative biomarker genes in gastric cancer with correlated expression and copy number, and the last objective was to characterize cancers based on DNA copy number amplifications. DNA microarrays, an in vitro model and real-time polymerase chain reaction were used to measure gene expression changes in H. pylori infected AGS cells. In order to identify the transcription factors and signal transduction pathways that were activated after H. pylori infection, gene expression profiling data from the H. pylori experiments and a bioinformatics approach accompanied by experimental validation were used. Genome-wide expression and copy number microarray analysis of clinical gastric cancer samples and immunohistochemistry on tissue microarray were used to identify putative gastric cancer genes. Data mining and machine learning techniques were applied to study amplifications in a cross-section of cancers. FOS and various stress response genes were regulated by H. pylori infection. H. pylori regulated genes were enriched in the chromosomal regions that are frequently changed in gastric cancer, suggesting that molecular pathways of gastric cancer and premalignant H. pylori infection that induces gastritis are interconnected. 16 transcription factors were identified as being associated with H. pylori infection induced changes in gene expression. NF-κB transcription factor and p50 and p65 subunits were verified using elecrophoretic mobility shift assays. ERBB2 and other genes located in 17q12- q21 were found to be up-regulated in association with copy number amplification in gastric cancer. Cancers with similar cell type and origin clustered together based on the genomic localization of the amplifications. Cancer genes and large genes were co-localized with amplified regions and fragile sites, telomeres, centromeres and light chromosome bands were enriched at the amplification boundaries. H. pylori activated transcription factors and signal transduction pathways function in cellular mechanisms that might be capable of promoting carcinogenesis of the stomach. Intestinal and diffuse type gastric cancers showed distinct molecular genetic profiles. Integration of gene expression and copy number microarray data allowed the identification of genes that might be involved in gastric carcinogenesis and have clinical relevance. Gene amplifications were demonstrated to be non-random genomic instabilities. Cell lineage, properties of precursor stem cells, tissue microenvironment and genomic map localization of specific oncogenes define the site specificity of DNA amplifications, whereas labile genomic features define the structures of amplicons. These conclusions suggest that the definition of genomic changes in cancer is based on the interplay between the cancer cell and the tumor microenvironment.

Genomic alterations in Myeloproliferative Neoplasms and Myelodysplastic Syndromes

Myeloproliferative neoplasms (MPN) and myelodysplastic syndromes (MDS) are a heterogeneous group of clonal hematopoietic disorders whose etiology and molecular pathogenesis are poorly understood. During the past decade, enormous developments in microarray technology and bioinformatics methods have made it possible to mine novel molecular alterations in a large number of malignancies, including MPN and MDS, which has facilitated the detection of new prognostic, predictive and therapeutic biomarkers for disease stratification. By applying novel microarray techniques, we profiled copy number alterations and microRNA (miRNA) expression changes in bone marrow aspirate and blood samples. In addition, we set up and validated an miRNA expression test for bone marrow core biopsies in order to utilize the large archive material available in many laboratories. We also tested JAK2 mutation status and compare it with the in vitro growth pattern of hematologic progenitors cells. In the study focusing on 100 MPN cases, we detected a Janus kinase 2 (JAK2) mutation in 71 cases. We observed spontaneous erythroid colony growth in all mutation-positive cases in addition to nine mutation negative cases. Interestingly, seven JAK2V167F negative ET cases showed spontaneous megakaryocyte colony formation, one case of which also harbored a myeloproliferative leukemia virus oncogene (MPL) mutation. We studied copy number alterations in 35 MPN and 37 MDS cases by using oligonucleotide-based array comparative hybridization (array CGH). Only one essential thrombocythemia (ET) case presented copy number alterations in chromosomes 1q and 13q. In contrast, MDS cases were characterized by numerous novel cryptic chromosomal aberrations with the most common copy number losses at 5q21.3q33.1 and 7q22.1q33, while the most common copy number gain was trisomy 8. As for the study of the bone marrow core biopsy samples, we showed that even though these samples were embedded in paraffin and underwent decalcification, they were reliable sources of miRNA and suitable for array expression analysis. Further, when studying the miRNA expression profiles of the 19 MDS cases, we found that, compared to controls, two miRNAs (one human Epstein-Barr virus (miR-BART13) miRNA and one human (has-miR-671-5p) miRNA) were downregulated, whereas two other miRNAs (hsa-miR-720 and hsa-miR-21) were upregulated. However, we could find no correlation between copy number alterations and microRNA expression when integrating these two data. This thesis brings to light new information about genomic changes implicated in the development of MPN and MDS, and also underlines the power of applying genome-wide array screening techniques in neoplasias. Rapid advances in molecular techniques and the integration of different genomic data will enable the discovery of the biological contexts of many complex disorders, including myeloid neoplasias.

Characterization of a genomic signature of pregnancy identified in the breast.

The objective of this study was to comprehensively compare the genomic profiles in the breast of parous and nulliparous postmenopausal women to identify genes that permanently change their expression following pregnancy. The study was designed as a two-phase approach. In the discovery phase, we compared breast genomic profiles of 37 parous with 18 nulliparous postmenopausal women. In the validation phase, confirmation of the genomic patterns observed in the discovery phase was sought in an independent set of 30 parous and 22 nulliparous postmenopausal women. RNA was hybridized to Affymetrix HG_U133 Plus 2.0 oligonucleotide arrays containing probes to 54,675 transcripts, scanned and the images analyzed using Affymetrix GCOS software. Surrogate variable analysis, logistic regression, and significance analysis of microarrays were used to identify statistically significant differences in expression of genes. The false discovery rate (FDR) approach was used to control for multiple comparisons. We found that 208 genes (305 probe sets) were differentially expressed between parous and nulliparous women in both discovery and validation phases of the study at an FDR of 10% and with at least a 1.25-fold change. These genes are involved in regulation of transcription, centrosome organization, RNA splicing, cell-cycle control, adhesion, and differentiation. The results provide initial evidence that full-term pregnancy induces long-term genomic changes in the breast. The genomic signature of pregnancy could be used as an intermediate marker to assess potential chemopreventive interventions with hormones mimicking the effects of pregnancy for prevention of breast cancer.

Genomic signatures of evolutionary transitions from solitary to group living

The evolution of eusociality is one of the major transitions in evolution, but the underlying genomic changes are unknown. We compared the genomes of 10 bee species that vary in social complexity, representing multiple independent transitions in social evolution, and report three major findings. First, many important genes show evidence of neutral evolution as a consequence of relaxed selection with increasing social complexity. Second, there is no single road map to eusociality; independent evolutionary transitions in sociality have independent genetic underpinnings. Third, though clearly independent in detail, these transitions do have similar general features, including an increase in constrained protein evolution accompanied by increases in the potential for gene regulation and decreases in diversity and abundance of transposable elements. Eusociality may arise through different mechanisms each time, but would likely always involve an increase in the complexity of gene networks.

Genomic alterations in patients showing multiple primary tumors and family history of cancer

Multiple primary tumors (MPT) are a major cause of mortality and morbidity among patients that have survived after the treatment of a first cancer. It has been proposed that after the first primary tumor, high risk of a subsequent tumor could be associated with radiotherapy used as treatment for the first cancer. Other potential risk factors include unhealthy lifestyle, genetic predisposition, aging, environmental determinants or an interaction between these factors. However, an association between the presence of MPT and family history of cancer in cases without clinical and molecular evidence of a known hereditary cancer syndrome is rarely described. Genomic DNA from 12 patients with at least two primary tumors and without mutations on TP53 was evaluated by CytoScan HD Array (Affymetrix). Chromosome Analysis Suite (ChAS) software v.2.0.1 was used considering at least 50 markers for gains; 25 for losses and a minimum of 5Mb for cnLOHs. Data from 1038 phenotypically healthy individuals (Affymetrix) and from Database of Genomic Variants were used as reference. Only alterations found in <1% (rare) or never described (new rare) in the reference population were considered. All cases, except one, presented a family history of cancer. Five cases developed MTP after radiotherapy and only one was located in the same treated area. It was detected 67 rare and 15 new rare genomic alterations encompassing 5.906 genes: 17 losses, 29 gains, and 36 cnLOH. X chromosome presented the higher number of alterations. Two patients with breast cancer presented a large deletion/cnLOH on 7q21. Enrichment analysis revealed 1275 genes associated with breast cancer (p= 0.001), which was diagnosed in 6 patients and their family members (all negative for BRCA1/2 or TP53 mutations). cnLOHs accounted for 44% of all the alterations. A significant proportion of cases (11/12) presented family history of cancer and the patients were not submitted to radiotherapy (7/12). We demonstrated the presence of rare genomic alterations in patients with MPT suggesting their involvement in the MPT development. cnLOH may arise as a new mechanism associated with the risk to develop MPT. All authors have declared no conflicts of interest.

Genomic gains in prostatic carcinoma and proliferative inflammatory atrophy in dogs

The dog can spontaneously develop prostate cancer and consequently can be used as an experimental model for prostatic diseases associated with aging, including benign prostate hyperplasia (BPH) and prostate carcinoma (PCa). DNA copy number variations (CNVs) have been used to identify genes associated with cancer development and progression. DNA microarray based comparative genomic hybridization (aCGH) is a technique that allows to identify copy number of thousands of genes throughout the genome. aCGH was used to identify genomic regions with significantly different DNA copy number in three benign prostatic hyperplasia (BPH), four proliferative inflammatory atrophy (PIA), and 14 canine prostate carcinoma (PCa). Five histologically normal prostate tissue were used as reference. Genomic DNA was extracted from formalin fixed and paraffin embedded samples and CNVs data was evaluated in Canine Genome CGH Microarray 4x44K (G2519F, Design ID021193, Agilent). Data analysis was performed using Genomic Workbench Standard Edition 5.0.14 (Agilent). PCa showed higher number of altered genes related to canonical diseases process, cellular functions and molecular pathways as well as greater inter-relationship between genes, compared with PIA and BPH. In conclusion, PCa showed a more complex genotype, being losses the most frequent genomic changes. Some discrepancies between genomic alterations in human and canine carcinomas may indicate the different clinical behavior of these tumors in these two species. In addition, it was observed was an ascending pattern of genomic complexity in BPH, PIA and CA consistent with a model of multistep tumor progression.

Microarray-based comparative genomic hybridisation of breast cancer patients receiving neoadjuvant chemotherapy

We analysed the molecular genetic profiles of breast cancer samples before and after neoadjuvant chemotherapy with combination doxorubicin and cyclophosphamide (AC). DNA was obtained from microdissected frozen breast core biopsies from 44 patients before chemotherapy. Additional samples were obtained before the second course of chemotherapy (D21) and after the completion of the treatment (surgical specimens) in 17 and 21 patients, respectively. Microarray-based comparative genome hybridisation was performed using a platform containing approx5800 bacterial artificial chromosome clones (genome-wide resolution: 0.9 Mb). Analysis of the 44 pretreatment biopsies revealed that losses of 4p, 4q, 5q, 12q13.11–12q13.12, 17p11.2 and 17q11.2; and gains of 1p, 2p, 7q, 9p, 11q, 19p and 19q were significantly associated with oestrogen receptor negativity. 16q21–q22.1 losses were associated with lobular and 8q24 gains with ductal types. Losses of 5q33.3–q4 and 18p11.31 and gains of 6p25.1–p25.2 and Xp11.4 were associated with HER2 amplification. No correlations between DNA copy number changes and clinical response to AC were found. Microarray-based comparative genome hybridisation analysis of matched pretreatment and D21 biopsies failed to identify statistically significant differences, whereas a comparison between matched pretreatment and surgical samples revealed a statistically significant acquired copy number gain on 11p15.2–11p15.5. The modest chemotherapy-driven genomic changes, despite profound loss of cell numbers, suggest that there is little therapeutic selection of resistant non-modal cell lineages.